Experimental test for subdominant superconducting phases with complex order 
parameters in cuprate grain boundary junctions 
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We propose and implement a direct experimental test for subdominant superconducting phases 
with broken time-reversal symmetry in d-wave superconductors. The critical current of 45°- 
asymmetric grain boundary junctions are shown to be extremely sensitive to the predicted onset 
of a complex order parameter at (llO)-surfaces and near magnetic impurities. Measurements in 
YBCO and Ni-doped YBCO junctions indicate that the symmetry at the surface is consistent with 
pure d-wave at all temperatures, putting limits on the magnitude and chiral domain structure of 
any subdominant symmetry component. 

PACS numbers: 74.20.Rp, 74.50.+r 



It has recently been recognized that the order parame- 
ters of unconventional superconductors may be unstable 
in the presence of perturbations as a result of their strong 
magnitude and phase anisotropy. In the high tempera- 
ture cuprate superconductors, the d x 2_ y 2 state dominant 
in the bulk |j],E| is readily suppressed at surfaces, at in- 
terfaces with other materials, in vortex cores, and in the 
vicinity of impurities. This suppression may allow the 
emergence of localized regions with different symmetries, 
possibly including phases with complex order parameters 
that break time-reversal symmetry. These states are of 
scientific interest, offering unique opportunities for study- 
ing novel phases in superconductor systems, and may also 
be important for understanding the microscopic pairing 
mechanism and for the implementation of high tempera- 
ture superconducting materials in electronic devices. 

The fragility of the d-wave state is a result of the 
phase change of ir between orthogonal directions. At 
a (llO)-surface, all specular reflection trajectories un- 
dergo a phase change of tt, forming via Andreev reflection 
bound surface states with zero energy ||. These states 
arc observed as a zero bias conductance peak in tunnel- 
ing spectroscopy experiments and as a modifica- 
tion of the low temperature penetration depth j?]]^] . The 
occupation of these states suppresses the d-wave order 
parameter, and it has been suggested that a secondary 
pairing interaction or subdominant component of the pri- 
mary pairing interaction could induce a complex mixture 
phase at the surface |^|[(|- The observation of a zero- 
magnetic field splitting of the zero bias conductance peak 
in a-b plane quasiparticle tunneling [^,[nj, may be evi- 
dence for broken time-reversal symmetry, consistent with 
a complex superconducting order parameter at the sur- 
face. A similar effect may occur due to scattering from 
magnetic impurities in the cuprates [O], and some ex- 
periments have observed an abrupt drop in the thermal 
conductivity of Ni-doped BSCCO crystals M] that may 
be interpreted as the opening of an energy gap in all di- 



rections, again consistent with the onset of a complex 
order parameter. 

In this Letter, we present an experiment designed to 
test specifically for the formation of a complex order 
parameter along the (llO)-surface of YBCO films by 
measuring the variation of the critical current of grain 
boundary junctions with temperature and magnetic field. 
We present calculations demonstrating that the onset of 
a complex superconducting order parameter has a dra- 
matic effect on the magnitude and magnetic field diffrac- 
tion patterns of the supercurrent in this geometry. Mea- 
surements of the critical current in 45°-asymmetric grain 
boundary junctions of YBCO over a wide temperature 
range are consistent with pure d x 2_ y 2 symmetry, show- 
ing no evidence for complex secondary phases. Similar 
results are obtained for Ni-doped films in which a bulk 
transition to a complex state has been suggested. 

The most straightforward way to test for the pres- 
ence of a complex order parameter is to perform a cor- 
ner SQUID or corner junction experiment with planar 
faces orthogonal to the (100) and (llO)-directions [0. 
In this configuration, a complex order parameter has a 
phase shift between and tt, yielding a characteristic 
modulation pattern. However, since smooth (110) faces 
are not attainable in cuprate crystals, we instead make 
use of the extensively-studied bicrystal thin film grain 
boundary Josephson junction. The critical current of 
such junctions is a strong function of the misorientation 
angle 9, falling off exponentially fTi) to a minimum value 
at 8 = 45° that is roughly three orders of magnitude 
smaller than the thin film critical current density. This 
variation is largely due to the due to the anisotropy of 
the d-wave order parameter, with additional contribu- 
tions from the structural mismatch at the interface and 
band-bending effects Another important feature of 
grain boundary junctions is faceting of the barrier plane 
caused during film growth by the competing grain ori- 
entations along the interface. For most grain boundary 
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orientations, the faceting has little effect and the junc- 
tion critical current varies with magnetic field according 
to the usual Fraunhofer diffraction pattern. However, for 
45°-asymmetric junctions, in which the lobe of the bulk 
d-wave order parameter in one electrode and the node 
in the other arc oriented normal to the interface, the 
critical current modulates with applied magnetic field in 
a complicated manner that reflects the detailed faceting 
along the grain boundary and the symmetry of the order 
parameter at the junction interface |l6[ |. 

To demonstrate this, we calculate the critical current 
of the grain boundary junction as a function of applied 
magnetic field for different order parameter symmetries. 
To model the grain boundary junction, we assume a ran- 
dom distribution of facet widths and angles. This distri- 
bution is unique to each junction, depending strongly on 
the angle and uniformity of the substrate grain boundary 
interface and the growth kinetics of the deposited cuprate 
film. From AFM images, it is estimated that typical facet 
widths range from 5 — lOOnm, with a wide range of an- 
gles. We assume that the local critical current density 
is proportional to the product of the magnitudes of the 
order parameters in the two electrodes and the relative 
phase across the junction, which includes the local dif- 
ference between the quantum mechanical phases in the 
electrodes resulting from the unconventional symmetry 
as well as that produced by the magnetic flux threading 
the barrier. Because of faceting, these quantities depend 
on the local orientation of the grain boundary interface, 
which selects the tunneling direction, relative to the a- 
b axes of the films in the electrodes. We consider only 
the effect of the applied magnetic field, neglecting any 
fields produced by the tunneling currents. This "short 
junction" limit is valid since the width of our junctions 
(5 — 20/im) is small compared to the Josephson penetra- 
tion depth A j = ($o/2ttij,otJ c ) 1 / 2 , typically 25 — bOfim 
for our devices. 

Shown in Figure 1 is the meander profile for a 10/im- 
wide junction with 25 facets randomly-angled from +45° 
to —45° from the nominal grain boundary interface. We 
calculate the critical current for this junction as a func- 
tion of the magnetic flux threading the barrier for four 
pairing symmetries: s, d x 2_ y 2 1 and the complex mixture 
states d x 2_ y 2 + id xy and d x 2_ y 2 + is. For s-wave, the 
magnitude and phase of the order parameter is uniform, 
giving the Fraunhofer diffraction modulation pattern fa- 
miliar from single slit optical interference and typically 
observed for Josephson junctions in conventional BCS 
superconductors. For d-wave symmetry, the variation 
of facet angles causes a corresponding variation of the 
magnitude and, more importantly, the sign of the or- 
der parameter across the junction, resulting in a sub- 
stantial suppression of the critical current and a com- 
plicated modulation with applied field that persists out 
to high magnetic fields. There are several characteris- 
tic features of the field modulation pattern for d-wave 



symmetry. First, the maximum critical current always 
occurs at a finite magnetic field for which the various 
facets of the junction are brought into phase. Second, 
despite its complex structure, the modulation pattern is 
always symmetric with respect to magnetic field polarity 
in the short junction limit; including self-field effects can 
introduce a small polarity asymmetry due to the inho- 
mogeneous current flow across the junction. Finally, the 
fundamental field modulation period observed is set by 
the overall width and magnetic barrier of the junction 
and so is the same for all random facet distributions, 
whereas the critical current at low fields may exhibit ei- 
ther a peak or a dip, depending on the faceting. 
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FIG. 1. (a) Faceting for the 45°-asymmetric grain bound- 
ary interface, (b) Calculated critical current modulation for 
this junction for s, d x 2_ y 2, d x 2_ y 2 + id xy and d x 2_ y 2 + is 
symmetries, assuming e = 0.20. 

For a complex order parameter, d + id' or d + is, the 
modulation pattern is significantly altered. The zero field 
critical current is dramatically enhanced, which is ex- 
pected since the presence of the secondary pairing com- 
ponent removes the node in the order parameter facing 
the interface in one of the electrodes, allowing a substan- 
tial tunneling contribution. More strikingly, the polarity 
symmetry is broken, resulting in an asymmetry in the 
critical current modulation pattern that is is most pro- 
nounced for small magnetic fields. This asymmetry is 
a direct manifestation of the broken time-reversal sym- 
metry characteristic of a complex superconducting order 
parameter. For different facet distributions, the detailed 
shape of the diffraction pattern changes, but the key qual- 
itative features remain the same. 
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These phenomena provide a sensitive test for the exis- 
tence of a complex order parameter in the superconduct- 
ing electrodes of a grain boundary junction. In Figure 
2(a), we demonstrate the effect of an s-wave secondary 
out-of-phase order parameter component by showing crit- 
ical current diffraction patterns for several values of e, 
the relative size of the subdominant s component in a 
d x 2_ y 2 + is superconductor. To quantify the sensitiv- 
ity, Figure 2(b,c) shows the enhancement of the zero-field 
critical current I c (0) and the onset of the fractional polar- 
ity asymmetry a(H) = {[I c {+H) - I c (-H)]/[I e (+H) + 
I c (—H)]}, averaged over values of applied flux from 
— IO$o to +10<I>o, as a function of e for d + is and d + id' 
symmetries. If a secondary order parameter component 
onsets at a specific phase transition temperature well be- 
low T c , as has been predicted, we would expect e to on- 
set and grow as the temperature is decreased. Figure 
2(d) shows the critical current vs. temperature predicted 
for the modeled junction, using a simple model in which 
I C (T) ~ A(T), which we assume to have a BCS tempera- 
ture dependence. This sharp increase in critical current, 
and the simultaneous onset of polarity asymmetry, would 
be a definitive signature of the nucleation of a complex 
order parameter. 
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FIG. 2. (a) Calculated critical current modulation of the 
junction in Figure 1 for d x 2_ y 2 + is symmetry and a range 
of subdominant components < e < 1. (b) Calculated en- 
hancement of -Zc(O) vs. e for d x 2_ y 2 + id xy and d x a_ y a + is 
symmetries, (c) Calculated onset of critical current asymme- 
try a{H) = {[Ic(+H) - I C (-H)]/[I C (+H) + /c(--ff)]} vs. e 
for d x 2_ y 2 +id xy and d x 2_ y 2 +is symmetries, (d) Simulated 



critical current vs. temperature for s, d x 
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+ is symmetries, assuming e = 0.20. 



We have made measurements of the magnetic field vari- 
ation of the critical current of a large number of 45°- 
asymmetric grain boundary junctions over a wide range 
of temperatures from T c down to 0.3K. YBa2Cu307_ x 
thin films of thickness lOOnm were grown on 45°- 
asymmetric bicrystal substrates by pulsed laser deposi- 
tion (growth temperature 830C, oxygen pressure 0.5Torr, 
power 450mJ), yielding a T c of 85-92K, with a transi- 
tion width of 1-2K, as determined from two-coil mag- 
netic screening measurements. We also made films with 
Ni doping (3% -5%), which had a slightly lower T c ( 80K) 
and a broader transition ( 5K). The films were patterned 
into strips of width 5 — 20/im using Ar-ion milling, pro- 
ducing grain boundary junctions with typical areas of 
I0 _8 cm 2 . The Ni-doped films were measured in a dilu- 
tion refrigerator down to lOOmK. 

To measure the critical current, we used a feedback 
technique in which the bias current is automatically con- 
trolled to maintain a small voltage level across the junc- 
tion, typically 10 [iV, as the magnetic field or temperature 
is varied. The critical currents increase monotonically as 
the temperature is lowered below T c , exhibiting a nearly 
linear dependence over much of the temperature range, 
as shown in the inset of Figure 3. The zero field critical 
current is typically 1 — 10fiA, corresponding to an aver- 
age current density of 100 — 1000 A/ cm 2 . The modulation 
of the critical current at T=4K for this junction, plotted 
in Figure 3, has the complex structure expected for a 
junction with pure d-wave symmetry in which the order 
parameter changes sign. The critical current is maximum 
at about 12G, rather than at zero magnetic field, and is 
nearly symmetric with respect to polarity. The measured 
average asymmetry parameter (a) at the peak currents is 
no more than 3% for this junction, typical of all junctions 
we have measured. 
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Measured critical current modulation at T — 4K 
for a 45°-asymmetric YBCO grain boundary interface. Inset 
shows the measured / c (0) vs. temperature . 



FIG. 3. 



Figure 4 shows the critical current vs. field at a series 
of different temperatures for (a) a pure YBCO, and (b) 
a 3%-Ni-doped junction. In both cases, the magnitude 
of the critical current increases with decreasing tempera- 
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ture, but the modulation pattern retains the same shape, 
reproducing each detailed feature, and in particular re- 
mains symmetric with respect to field polarity over the 
entire temperature range. 
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FIG. 4. Measured critical current modulation at a series 
of temperatures for: (a) a pure YBCO junction, and (b) a 
5%-Ni doped YBCO junction. 

In all junctions studied to date, both pure and Ni- 
doped, the polarity asymmetry of junctions cooled in 
zero field is typically a few percent and is never larger 
than 10%. This level of asymmetry is consistent with 
that expected from experimental uncertainty, self-field ef- 
fects, and/or trapped magnetic flux in the vicinity of the 
junction. More significantly, we have never observed any 
abrupt increases in either the zero field critical current 
or the polarity asymmetry, which have been predicted 
to occur at a phase transition to a state with a com- 
plex order parameter as the junctions are cooled. Based 
on a comparison of our simulations and measurements, 
we would place an upper bound on the fractional mag- 
nitude of an out-of-phase s or d xy component added to 
the d x 2_ y 2 order parameter at about 1%. We note that 
we also do not find evidence for the onset of a real (in- 
phase) subdominant order parameter component, which 
would exhibit an increased zero-field critical current but 
no polarity asymmetry. 

It is important to consider possible reasons why we 
might not observe the onset of a substantial complex or- 
der parameter, as has been predicted by theoretical treat- 
ments. One scenario is that the complex order parameter 
region induced at the interface forms domains of alter- 
nating chirality, e.g. d + is or d — is. Such domaining 



is not energetically favorable due the domain wall en- 
ergy (of order the Josephson coupling energy), but could 
nucleate as metastable states when a secondary phase 
onsets. We have simulated the effect of chirality domain- 
ing by putting a random spatial distribution of domains 
along the interface on the (llO)-electrode of the junc- 
tion. Because of the alternating sign of the secondary 
order parameter component, the net increase in the criti- 
cal current in the complex phase is substantially reduced. 
However, the domaining effectively increases the density 
of facets on which the order parameter phase modulates, 
increasing both the magnetic field range and the magni- 
tude of the polarity asymmetry. Thus, domain formation 
should enhance rather than obscure the formation of a 
complex order parameter at the interface. 

A second scenario is that no complex superconducting 
phases are formed. One reason could be that the rel- 
atively high transmission coefficient of grain boundary 
junctions does not allow sufficient Andreev reflection to 
produce the zero energy bound states that are responsi- 
ble for the suppression of the d-wave order parameter at 
the surface. However, it should be noted that the trans- 
mission coefficient of the 45°-asymmetric grain bound- 
ary junctions considered here is not particularly large, 
with current densities of only 100 — 1000A/cm 2 . This is 
supported by measurements of the conductance vs. volt- 
age in such junctions which show a pronounced zero bias 
conductance peak |l7| , similar to that observed in planar 
tunnel junctions on thin films and crystals with in-plane 
orientations. We also note that is has been predicted 
that proximity coupling between the electrodes may in- 
duce complex phases even in high transmission junctions 
p8| . A second option is that the d-wave order param- 
eter is suppressed, but that the conditions near the in- 
terface (density of states, pairing interaction, ...) are not 
conducive to the formation of a secondary out-of-phase 
superconducting component so that the d x i_ y i state re- 
mains stable at all temperatures. In either case, the re- 
sults suggest that the apparent time-reversal symmetry 
breaking observed in tunneling spectroscopy, low temper- 
ature transport, and local magnetic field measurements 
may arise from some microscopic mechanism, rather from 
than nucleation of a macroscopic superconducting region 
with complex order parameter symmetry. Primary candi- 
dates include magnetic surface states, antiferromagnetic 
bond currents, and barrier defects. 

In conclusion, we have proposed and implemented a 
new experimental test for the onset of complex order pa- 
rameter symmetry at surface and near impurities in un- 
conventional superconductors. The critical current mag- 
nitude and magnetic field modulation of 45°-asymmetric 
grain boundary junctions are shown to be extremely sen- 
sitive to the onset of an out-of-phase secondary order 
parameter, characterized by a dramatic enhancement in 
the zero field current and an asymmetry with respect to 
magnetic field polarity. Measurements in YBCO and Ni- 
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doped YBCO junctions are consistent with pure d x 2_ y 2 
symmetry, with less than 1% mixture of a subdominant 
superconducting phase. 

We wish to thank Joe Hilliard, Chris Michael, and 
Tony Banks for vital technical assistance. This work 
is supported by the National Science Foundation grant 
NSF-DMR99-72087 and by the Department of Energy 
grant DEFDG02-96ER45439. We acknowledge extensive 
use of the Microfabrication Laboratory and the DOE 
Center for Microanalysis of Materials in the Frederick 
Seitz Materials Research Laboratory. 



[1] D. A. Wollman et al., Phys. Rev. Lett. 71, 2134 (1993); 

D. J. Van Harlingen, Rev. Mod. Phys. 67, 515 (1995). 
[2] C. C Tsuei et al., Phys. Rev. Lett. 73, 593 (1994). 
[3] C. R. Hu, Phys. Rev. Lett. 72, 1526 (1994). 



[4] J. Geerk et al., Z. Phys. 73, 329 (1988). 
[5] J. Lesueur et al., Physica C 191, 325 (1992). 
[6] M. Covington et al., Phys. Rev. Lett. 79, 277 (1997). 
[7] H. Walter et al., Phy s. Rev. Lett. 80, 35 98 (1998). 
[8] A. Carrington et al., |cond-mat/0011182; . 
[9] M. Matsumoto and H. Shiba, J. Phys. Soc. Japan 64, 
1703 (1995). 

[10] L. J. Buchholtz et al., J. Low Temp. Phys. 101, 1079 
(1995). 

[11] R. Krupke and G. Deutscher, Phys. Rev. Lett. 83, 4634 
(1999). 

[12] A.V. Balatsky, Phys. Rev. Lett. 80, 1972 (1998). 
[13] R. Movshovich et al., Phys. Rev. Lett. 80, 1968 (1998). 
[14] D. Dimos et al., Phys. Rev. Lett. 61, 219 (1988). 
[15] H. Hilgenkamp and J. Mannhart, Appl. Phys. Lett. 73, 
265 (1998). 

[16] H. Hilgenkamp, J. Mannhart, and B. Mayer, Phys. Rev.B 

53, 14586 (1996). 
[17] L. Alffet al., Phys. Rev. B 55, 14757 (1997). 
[18] M. Fogelstrom and S.-K. Yip, Phys. Rev. B 57, R14060 

(1998). 



5 



